B-mesons from top-quark decay in presence of the charged-Higgs 
boson in the Zero-Mass Variable-Flavor-Number Scheme 
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We study the energy spectrum of the inclusive bottom-flavored mesons in top-quark decay into a 
charged-Higgs-boson and a massless bottom quark at next-to-leading order (NLO) in the two-Higgs- 
doublet model. To extract the resuh we work in the Zero-Mass Variable-Flavor-Number scheme(ZM- 
VFNs) using reahstic nonperturbative fragmentation functions obtained through a global fit to e'^ e~ 
data from CERN LEPl and SLAC SLC on the Z-boson resonance. We study both the contribution 
of the bottom-quark fragmentation and the gluon fragmentation to produce the bottom-fiavored 
meson (B-Meson). We find that the contribution of the gluon leads to an appreciable reduction in 
decay rate at low values of the B-meson energy. It means the NLO corrections are significant. 

PACS numbers: 12.38.Bx, 13.85.Ni, 14.40.Nd, 14.65.Ha, 14.80.Da 



I. INTRODUCTION 

Top quark is the heaviest elementary particle with a 
mass of 172.0 GeV [l|. Top's large mass is a reason to 
rapid decay so that it has no time to hadronize and if it 
were not for the confinement of color, the top quark could 
be considered as a free particle. This property allows the 
top quark to behave like a real particle and one can safely 
describe its decay in perturbative theory. The Large 
Hadron Collider (LHC) is a superlative top factory with 
90 million tt-pairs per year which will allow us to perform 
accurate studies of the top-quark properties, such as its 
mass mt, total decay width Tt, the effective weak mix- 
ing angle, and elements Vtq of the Cabibbo-Kobayashi- 
Maskawa (CKM) 01 quark mixing matrix. The theoreti- 
cal aspects of top-quark physics at the LHC are summa- 
rized in a recent paper [3|. 

Since \ Vtb\ ~ 1, top quarks almost exclusively decay to 
bottom quarks, via t — >■ hW^ within the Standard Model 
(SM) theory and in beyond-the-SM theories with an ex- 
tended Higgs sector, top quarks decay via t — )■ hH^ . 
Many extensions of the Standard Model contain more 
than one Higgs doublet, and the new degrees of freedom 
appear as extra Higgs scalars. For example, in the super- 
symmetric SM, at least two Higgs doublets are needed 
so as to cancel gauge anomalies and to generate masses 
for both up- and down- type quarks in the Weinberg 
model of CP violation, at least three Higgs doublets are 
needed in order to have spontaneous CP violation(see [5| 
and references therein). Both neutral and charged phys- 
ical Higgs bosons exist in all these extensions of the SM. 
The observation of charged Higgs bosons, H^, would in- 
dicate physics beyond the SM. For the sake of simplicity, 
we will limit ourselves throughout this paper to the decay 
t bH^ in a model with two-Higgs-doublet, in which 
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case there is only one physical charged Higgs boson re- 
maining after spontaneous symmetry breaking. New re- 
sults of a search for the charged Higgs bosons in proton- 
proton collision at a center-of-mass energy of \/s = 7 
TeV are reported by the ATLAS Collaboration [g where 
the T-|-jets channel in tt decays is used with a hadroni- 
cally decaying r lepton in the final state. The reported 
data leads to a limit on the product of branching ratios 
BR{t bH^) X BR{H^ Tv) of 0.03-0.10 for 
masses in the range 90 GeV< m/f± < 160 GeV. 
As it is mentioned in Ref . 0| , a clear separation between 
the t — > bW^ and t — > bH^ can be achieved in both 
the tlx pair production and the t/tX single top produc- 
tion at the LHC. The QCD corrections to the decay rate 
t ^ 6 + are known at next-to-next-to-leading order 
Q and the NLO electroweak corrections were found in 
Ref. 0. In this paper, we evaluate the first order QCD 
corrections to the decay of an unpolarized top quark into 
a charged Higgs boson. 

Since bottom quarks hadronize, via 6 — > B + X, there- 
fore the decay process t — >■ BH^ -t- X is of prime im- 
portance and the particular purpose of this paper is to 
make predictions for the energy spectrum of b-flavored 
mesons in top decay. This measurement will be impor- 
tant for future tests of the Higgs coupling in the mini- 
mal supersymmetric standard model (MSSM) at LHC. 
The hadronization of the bottom quark was considered 
in the NLO QCD analyses of the decay t — >■ hW^ in 
Refs. [l3,[n|. As it is shown in [l2|, the hadronization 
of the bottom quark is identified to be the largest source 
of uncertainty in measurement of the top-quark mass. 
To study the distributions in the B-meson scaled-energy 

, we employ the massless scheme or zero-mass variable- 
fiavor-number (ZM-VFN) scheme [l^ in the top-quark 
rest frame. In this scheme, the zero mass parton 
approximation is also applied to the bottom quark 
and the non-zero value of the b-quark mass only en- 
ter through the initial condition of the nonperturba- 
tive fragmentation function(non-FF). Nonperturbative 
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FF is describing the hadronization process b ^ B + X 
and is subject to Dokshitzer-Gribov-Lipatov-Alteralli- 
Parisi(DGLAP) evolution and it is scale dependent. 

This paper is organized as follows. In SecUll we study 
the inclusive production of a meson from top-quark decay 
considering the factorization theorem and DGLAP equa- 
tions. In Sec, mil we present the calculation of the 0{as) 
QCD corrections to the tree-level rate of t — > bH^. We 
work in ZM-VFN scheme neglecting the b-quark mass in 
our QCD corrections but will retain the arbitrary value 
of 171^+ ■ In Sec. IIVI we present our numerical analysis. 
In Sec. |Vl our conclusions are summarized. 



II. FORMALISM 

We study the inclusive production of a B-meson from 
top-quark decay 



t^b + H+{g) BH+ + X, 



(1) 



providing that the top-quark mass rrit, bottom-quark 
mass mil and the charged Higgs boson mass m'^ satisfy 
mt > mil + 771^ . The gluon in Eq. ([ij contributes to 
the real radiation at NLO and both the b quark and the 
gluon may hadronize to the B meson. In the equation 
above, X stands for the unobserved final state. 

If we denote the four-momenta of top-quark, b quark, 
gluon and B meson by Pt,Pb,Pg and pb, respectively, 
therefore in the top-quark rest frame the b quark, gluon, 
and B meson have energies Ei = pt ■ pi/mt{i = b, g, B), 
where ms < Eb < [mj + — m'^^)/{2mt) and 
< {El,, Eg) < {ml - m^j + )/{2mt). It is convenient 
to introduce the scaled energy fractions Xi = Ei/E™^^ 
{i = b,g,B). 

We wish to calculate the partial decay width of pro- 
cess (HI) differential in xb, dT/dxB, at NLO in the ZM- 
VFN scheme. According to the factorization theorem of 
the QCD-improved parton model [l^, the energy dis- 
tribution of a hadron B can be expressed as the convo- 
lution of the parton-level spectrum with the nonpertur- 
bative fragmentation function Da{z, fip), describing the 
hadronization a ^ B. 



dT 
dxB 



E 

a=b,g 



dXa dVa , f^B 



(2) 

where dVa/dxa is the parton-level differential width of 
the process t — > a-t- AT, with X comprising the boson 
and any other parton. Here, and /i^ are the factor- 
ization and the renormalization scales, respectively. At 
NLO, the scale /i/j is associated with the renormalization 
of the strong coupling constant. In principle, one can use 
two different values for the factorization and renormal- 
ization scales; however, a choice often made consists of 
setting jjLR = fJ-F and we shall adopt this convention for 
most of the result which we shall show. 
In next section, we present our analytic results for 
dta/dxa{a = b,g) at NLO in ZM-VFNS. 



III. ANALYTIC RESULTS FOR dtc,/dx^ 

A. Tree-Level Rate of t ^ bH+ in ZM-VFNS 

The coupling of the charged Higgs boson to the top 
and bottom quark in the minimal supersymmetric stan- 
dard model(MSSM) can either be expressed as a super- 
position of scalar and pseudoscaler coupling factors or as 
a superposition of right- and left-chiral coupling factors 
[l6| . Adopting the first approach, the Born term ampli- 
tude for the process t ^ b + iJ+ can be parametrized as 
A/o = Ub{a + b"f5)ut, and the second scheme leads to the 
Born amplitude Mq = Ub{gt{l + 75)/2 + .96(1 - j5)/2}ut 
where a = {gt + gb)/'^ and b — {gt ~ gb)/'^- In a model 
with two Higgs doublets and generic coupling to all the 
quarks, it is difficult to avoid flavor-changing neutral cur- 
rents. We, therefore, limit ourselves to models that nat- 
urally stop these problems by restricting the Higgs cou- 
pling. As it is explained in Ref. [l6j . the first possibility 
is to have the doublet Hi coupling to all bosons and the 
doublet H2 coupling to all the quarks (model /). This 
leads to the coupling factors 



gu. 



2\/2mw 



2V2: 



my/ 



Vtb{mt - mt) cot/3, 
-Vtb{mt +mb) cot p. 



(3) 



The second possibility is to have the H2 couple to the 
right-chiral up-type quarks {ur, cr, Ir), and the Hi cou- 
ple to the right-chiral down- type quarks (model //). This 
model leads to the coupling factors 
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mw 



Vtb{mt cot/3 -f TOh tan/3), 
Vt&(mt cot /3 - TOfc tan ^). (4) 



In equations above, tan /3 = V2IV1 is the ratio of the 
vacuum expectation values of the two electrically neutral 
components of the two Higgs doublets and the weak cou- 
pling factor g^ is related to the Fermi's constant Gf by 
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The total decay width of t H^b at LO is given by 



mt{o? -h 6^) ml 
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-(1 + 



Ai(i,I^,!^)^ + . 

mi m^ ' 



2mhmt 



a^-b^ 



m^B+ + ^ 



where A(a, &, c) = + b'^ + (? — 2{ab + be + ca) is the 
Kallen function. In the limit of vanishing b-quark mass, 
a = 6 in model /, therefore the tree level decay width 
simplifies to 
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Y cot^ /3. 



(5) 
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For model // one has 



mbTTit 



52 



TOj cot (3 + mT tan^ (3 



(6) 



when the left-chiral couphng term, proportional to 
nib tan /?, can become comparable to the right-chiral cou- 
pling term rrit cot /3, one cannot therefore naively set 
nib = in all expressions. For example, if we take 
nib = 4.90 GeV, = 172.0 GeV, mjy+ = 120 GeV and 
tan /? « 10 thus the second term in the curly brackets 
can become as large as C(6%) in model //. In this pa- 
per we adopt, with good approximation, the Born term 
presented in Eq. ([5]) in both models when mf, — > 0, more 
detail can be found in Ref. (T7| . 

In the following, we discuss the calculation of the 
0{as) QCD corrections to the tree- level decay rate of 
t ^ b + and we present the parton-level expressions 
for dT{t BH+ + X)/dxB at NLO in the ZM-VFN 
scheme. 



B. Virtual Corrections 

In the ZM-VFN scheme, where = is put from 
the beginning, all singularities including the soft- and 
coUinear gluon emission are regularized by dimensional 
regularization in D ~ A — 2e space-time dimensions to 
become single poles in e, which are subtracted at factor- 
ization scale and absorbed into the bare FFs according 
to the modified minimal-subtraction scheme (MS). In 
this scheme, nib only sets the initial scale /i"^' = 0{mb) 
of the DGLAP evolution. 

Adopting the on-shell mass-renormalization scheme, the 
virtual one-loop corrections to the tbH^-veitex includes 
both IR- and UV-singularities. Therefore, the contri- 
bution of virtual corrections into the differential decay 
width normalized to the Born width, reads 



where, the mass renormalization constant of the top 
quark reads 



Snit as{^iB) ^ , 3 
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and from Ref. [18| , for the wave function renormalization 
constants we have 
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In Eqs. (jlO|) and (jlip , em and ejjy represent infra- 
red(IR) and ultra-violet (UV) singularities which arise 
from the collinear- and the soft-gluon singularities, re- 
spectively. Therefore, the real part of the one-loop vertex 
corrections A; is given by 



A, 



, l + 2y 

y 



ln(l - y) - 2Lz2(y) 



2F-I-1 . 4 1 49 TT^^ 
£uv (Jr 8 12 
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where, F = 21n(l 



V) 



5/2, 



\n{Anfj.%/ni1) + je 
Cf = (A^c - l)/(2A^c) = 4/3 for AT^ = 3 quark colors, 
and Li2{x) = — Jq (dt/t) ln(l — t) is the Spence function. 
All UV-divergences are canceled after summing all vir- 
tual corrections up but the IR-singularities are remaining 
which are now shown by e. The virtual corrections to the 
differential decay rate is then given by 



1 dr^"' 
fo dxb 
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|M"'^|25(1-X6), 
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GFm1\Vtb\'^ coi^ P 1^/2, the scaled mass y 

Sp^JM^oMloop + 



where, 

is y = ni'j^+/nif and 

M^^^^AIo). Following Ref. |5], the renormalized ampli- 
tude of the virtual corrections can be written as 



1/2 



Mioop = "6 (Act + Ai){a + b-i5)ut, 



(8) 



where Act stands for the counter term and A; arises from 
the one-loop vertex correction. Since we neglect the b 
quark mass, thus the counter term of the vertex consists 
of the top-quark mass renormalization and the wave func- 
tion renormalizations as 



Ac 



1 
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C. Real Corrections 

As it is explained in Ref. [l^, to cancel the IR- 
singularities of the virtual corrections, one needs to in- 
clude the real gluon emission, namely, t H^bg. To 
calculate the contribution of the real corrections, we start 
form the definition of decay rate. As before, to regulate 
the IR- divergences we work in Z) = 4 — 2e dimensions, 
therefore from the definition of decay rate, one has 



dV 



real 



, 2(4-_D) 
H'F 

2nit 



\M-eal^^dPS{pt,Pb,Pg,PH^), (13) 



where, the Phase Space element reads 



jD-1 



dPS 
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jD-1 



Ph+ 



{2n)D-^2Eb {2n)D-^2EH+ {2n 
X {2n)^d^ {pt -pb- Ph+ - Pg) 



i-D-l 



2E„ 



(14) 
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For simplicity, we choose the top-quark rest-frame and 
to calculate the differential rate dV^^""^ / dxb normalized 
to the Born width, we fix the momentum of b-quark in 
Eq. . To get the correct finite term in the normalized 
differential decay rate, the Born width Fq will have to be 
evaluated in the dimensional regularization at O(e^), i.e. 
fo ^ fo{l-e(F + f/2) + e2(FV2-|-F/2 + f7/8-7rV4)}. 
Thus, the contribution of the real gluon emission into the 
normalized differential decay width is then given by 



1 dtl""^ 



OLs 

2^ 



Cp[5{l-x,)[^-\{F+\) + ^ + 



-lny-f2Lz2(f-y)- — + -] 



(1 - Xb) + 

+F] -f 2(f 4 



Ixl - 8xb 



i+21n., + ^^^,^ 
'lii(f-a-b) 



1 - Xfc 



(15) 



where the plus distributions are defined as usual. 



According to this scheme, in order to get the MS coef- 
ficient functions we shall have to subtract from Eq. (fT6|) 
the 0{as) term multiplying the characteristic MS con- 
stant (— 1/e + 7_E — ln47r), therefore, we obtain 



dV 



MS 



dxh 



= r, 



o|(5(l-a;6)-f 



2tt 



Pj^°gHxb)\n^ + CFAixb)\} 
/i 771 J J 



dx 
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(17) 



where Pi" -* are the time-like splitting functions at leading 
order UA 



1 + 3,,^ 
(1-^ + 2^^'"^^ 



D. Analytic Results for Partial Decay Rates 

The NLO expression for dVb/dxb is obtained by sum- 
ming the Born term, the virtual one-loop and the real 
contributions. 

Since, the B meson can be also produced from the frag- 
mentation of the emitted real gluon, therefore, we also 
need the differential decay rate dTg/dxg in the ZM-VFN 
scheme. To calculate the dVg/dxg, as before, we start 
form Eq. p3p by fixing the momentum of gluon, there- 
fore, there will be no soft singularity. 
Our results list here 



and the coefficient functions A{xb) and B{xg) read 
A{xb) = (5(l-xb)[2i^ln(l-y)-2^1n2/ 



B{Xg) 



-2Li2(y)+2Li2(l-y)-4 



1 + 

il-Xb) 
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1 + xl - 
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V I - Xb 
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dTb 

dxb 
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dXg 
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-Cf 



7i=;-ln47r) + Bi{xg) 



(16) 

where Ai{xb) and Biixg) are free of singularities and, in 
the following, their functional form will be shown. 
To subtract the coUinear singularities remaining in 
Eq. (jl6p . we apply the modified minimal subtraction 
{MS) scheme where the collinear singularities are ab- 
sorbed into the bare EEs. This renormalizes the EEs and 
generates the finite terms of the form ln(TOj//i|,) in 
the differential decay rates. 



-ln(l -(1 -,).,)] +^^^-^^—^[ 
7xg + — -6-{l- y){lQxl - 8xg + 12) + 



Xg{l-yf{Axl-4xg + 7) 



(18) 



In this work we identify fin ~ fip = nit, so that in 
Eq. pT|) the terms proportional to ln(mj//i|,) vanish. 
Integrating dVb/dxb of Eq. (fT7|) over Xb{0 < Xb < 1), we 
obtain the renormalized decay rate 



F = En 1- 



2tt 



2w 2 
-^ln2/+(5--)hi(l-y) + 

1 - y y 
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2Li2{y)-2Li2{l~y)- 



(19) 



This result is in agreement with Refs. [20|, |2l| and the 
corrected version of [5|. As it is seen from Eq. ([T^ . in 
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the limit /rrit — )■ 0(= y — > 0) the total decay rate is 
finite and proportional to Gpm^. At the opposite limit 
where m^+I'mt — > 1(= y — > 1), due to the ln(l — y) 
singularity setting the bottom-quark mass to be zero is 
no longer a valid approximation for the differential decay 
rates and our results must be improved considering a 
massive b-quark. 



IV. NUMERICAL ANALYSIS 

We are now in a position to present our phenomeno- 
logical results by performing a numerical analysis. In the 
MSSM, the mass of the charged Higgs is strongly corre- 
lated with the other Higgs boson masses. The charged 
Higgs boson mass in the MSSM is restricted at tree- 
level by mH+ > mw, Ref. This restriction does 
not hold for some regions of parameter space after in- 
cluding radiative corrections. Therefore, two key phe- 
nomenological parameters in the charged Higgs searches 
are the Higgs mass mjy+ and tan /3, which are model- 
dependent. Searches of the charged Higsses over a good 
part of the mu+ — tan j3 plane in the MSSM is a pro- 
gram that still has to be carried out and this belongs 
to the LHC experiments. We adopt from Ref. [l| the 
present limit > 79.3 GeV obtained from LEP. 
We also adopt from Ref. [l| the input parameter val- 
ues Gf = 1.16637 X 10-5 GeV-2, m* = 172.0 GeV, 
nib 4.90 GeV, and = 5.279 GeV. We evaluate 
ctT'^'im) at NLO in the MS scheme using Eq. (8) of 
Ref. [10| , with ?i j = 5 active quark flavors and the typical 

QCD scale = 231.0 MeV. As mentioned before, the 

MS ... . . 

b-quark mass only enter through the initial condition of 
the nonperturbative fragmentation function. We employ 
the nonperturbative _B-hadron FFs that were determined 
at NLO in the ZM-VFN scheme through a joint fit [13 to 
e+e--annihilation data taken by ALEPH [ll] and OPAL 
[Hi at CERN LEPl and by SLD [2^ at SLAG SLC. 
Specifically, the power ansatz Df,(z,/i'^') = A^z"(l — zY 
was used as the initial condition for the b B FF at 
Idp^ = 4.5 GeV; while the gluon and light-quark FFs 
were generated via the DGLAP evolution. The fit yielded 
N = 4684.1, a = 16.87, and /3 = 2.628. 

To study the scaled-energy (xb) distribution of the 
bottom-flavored hadrons produced in top-quark decay, 
we consider the quantity dr{t —>■ BH^ + X)/dxB- In 
Fig. [1] we show our prediction for the size of the NLO 
corrections, by comparing the LO (dotted line) and NLO 
(solid line) results, and the relative importance of the 
b ^ B (dashed line) and g ^ B (dot-dashed line) 
fragmentation channels at NLO, taking tan/3 = 10 and 
mH+ = 120 GeV. The same NLO FFs are used for the LO 
result. Fig. [T] shows that the NLO corrections lead to a 
significant enhancement of the decay rate in the peak re- 
gion and above. Furthermore, the peak position is shifted 
towards higher values oi xb- The gluon fragmentation 
leads to an appreciable reduction in decay rate at low- 
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Figure 1: dV{t — >■ BH'^ -f X)/dxB as a function of xb in 
the ZM-VFN (mt = 0) scheme with = 120 GeV and 

tan /3 — 10. The NLO result (solid line) is compared to the 
LO one (dotted line) and broken up into the contributions 
due to 6 — >■ B (dashed line) and g ^ B (dot-dashed hue) 
fragmentation. 



Xb region, for xb ^ 0.3. For example, the gluon frag- 
mentation decreases the size of decay rate up to 43% at 
Xb = 0.12. For higher values of xb, the b ^ B contri- 
bution is dominant. As we explained in section |lTl the 
mass of B-meson is responsible for the appearance of the 
threshold at xb = "^.ttib /(mt(l — y)) = 0.12. 

In Fig. (2b) of Ref. [l7| , the unpolarized rate is plotted 
as a function of tan /3 for a sample value of mjj+ = 120 
GeV. It can be seen that the rate is quite small for tan /3 
values exceeding tan/3 = 2. 

In Fig. [21 we study the energy distribution of the B- 
meson in different values of the tan/3, i.e. tan /? = 1, 5, 10 
and 15, for =120 GeV. It can be seen that when 

tan (3 is increased the decay rate is decreased, as Fq is 
proportional to cot^ /3. 

In Fig. (2a) of Ref. [l3|, the unpolarized rate is also 
shown as a function oi y = rriH+fnit for tan/3 = 10. 
The functional behavior of the rate shows that the rate 
is largest when — > and it drops to zero when 
171^+ — )■ mt- 

Adopting the limit > 79.3 GeV from Ref. [l|, 

in Fig. [3] we study the energy distribution of the B- 
meson in different values of the Higgs boson mass, i.e. 
mH+ = 80, 100, 120 and 150 GeV, by fixing tan/^ = 10. 
As mentioned, the mass of B-meson creates the thresh- 
olds at xb = 0.08 (for mH+ = 80 GeV), xb = 0.09 (for 
mH+ = 100 GeV), xb = 0.12 (for 771h+ = 120 GeV) and 
XB^ 0.26 (for mH+ = 150 GeV). 

In comparison with the Born rate ~ 1.364 

GeV, the rate into a charged Higgs is generally quite small 
except for small tan/3 values. One finds equality of the 
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Figure 2: dr{t — >■ BW^ + X)/dxB as a function of a-s in 
different values of tan/3 = 1,5, 10 and 15, with mjj+ = 120 
GeV. When the values of tan /3 exceed tan (3 = 2, the decay 
rate becomes quite small [l7l |. 




Figure 3: dT{t — BH^ + X)/dxB as a function of xb in the 
ZM-VFN scheme with different values of = 80, 100, 120 

and 150 GeV (tan/3 = 10). Thresholds at xb are shown. 



rates into a and iJ+ only at tan /? = 0.56 for 171^+ — 
120 GeV. However, such a small tan (3 value is excluded 
by the indirect limits in the (m^± , tan (3) plane [l^ . 

In Fig. m we compared the energy distribution of the 
B-meson produced in decay modes t — > BW^ + X and 
t BH+ + X with tan/3 = 1, mH+ = 120 GeV, 
niiY = 80.399 GeV and the parton-level differential rates 
dtbit bW+)/dxb and dfg(t bW+)/dXg extracted 
from Ref. (l^. The thresholds appear ai xb = 0.08 (for 
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Figure 4: xb spectrum in top decay considering the decay 
modes t BW+ + X (solid line) and t BH+ + X (dashed 
line), taking 771^1,-+ = 80.399 GeV, = 120 GeV and 

tan/3 = 1. 



t BW+ + X) and xb = 0.12 ( for t BH+ + X). 
The total top-quark decay width is obtained by sum- 
ming the two partial widths r(^h_|.^y+ and Tf.^i,_^_B+ or- 
der by order. However, as Figs. ^ and ^ show the ra- 
diative corrections coming from the supersymmetric sec- 
tor depend extremely on the Higgs mass and tan /3, but 
Fig. [4] shows that the contribution of the decay mode 
t -I- 6 is always larger than the one coming from 

t iJ + + b, see more detail in Ref. ■ 

V. CONCLUSIONS 

Top-quark decays within the Standard Model are com- 
pletely dominated by the mode t — > + b due to 
\Vtb\ = 1 to a very high accuracy. In beyond-the-SM 
theories with an extended Higgs sector, top-quark de- 
cays can also be done via t + b. This charged 
Higgs boson has been searched for in high energy exper- 
iments, in particular, at LEP and the Tevatron but it 
has not been seen so far. To study the Higgs boson and 
new physics beyond the SM at LHC, as a superlative top 
factory, we need to understand the energy distribution of 
B-meson produced from top-quark decay. The domiirant 
decay mode t — )• BW'^ -f X was studied in Refs. [l^ and 
fist - In this paper we studied the scaled-energy {xb) dis- 
tribution of B-meson in i — BH'^ -I- X at NLO in the 
ZM-VFN scheme. We employed the nonperturbative B- 
meson FFs determined by a global fit [l^ of experimental 
data from Z factories, relying on their universality and 
scaling violations (27| . 

Comparison of future measurements of dT/dxB at the 
LHC with our NLO predictions will be important for fu- 
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ture tests of the Higgs coupling in the minimal super- 
symmetric SM(MSSM) and it will be the primary source 
of information on the B-meson FFs. 

Acknowledgments 

I would like to thank Professor Gustav Kramer for 
reading and improving the manuscript and also for im- 



portant discussions and comments. I would also like to 
thank Professor Bernd A. Kniehl for his helpful advices 
and oppening the doors of perturbative QCD to me. This 
work was supported by Yazd university and the Institute 
for Research in Fundamental Science (IPM). 



[1] K. Nakamura et al. (Particle Data Group), J. Phys. G 

37, 075021 (2010). 
[2] N. Cabibbo, Phys. Rev. Lett. 10, 531 (1963); [13 

M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 49, 

652 (1973). 

[3] W. Bernreuther, J. Phys. G 35, 083001 (2008). [14 

[4| J. F. Gunion and H. E. Haber, Nucl. Phys. B 272 (1986) 
1 [Erratum-ibid. B 402 (1993) 567[. 

[5] C. S. Li and T. C. Yuan, Phys. Rev. D 42, 
3088 (1990) [Erratum-ibid. D 47, 2156 (1993 
PHRVA,D47,2156.1993)]. [15 

[6[ ATLAS Collaboration, Search for charged Higgs bosons [16 
in the r+jets final state in tt decays with 1.03/6^^ of pp 
collision data recorded at y/s = 7 TeV with the ATLAS 
experiment, ATLAS-CONF-2011-138 (2011) [17 

[7[ A. AH, F. Barreiro and J. Llorente. .arXiv7ll03. 1827l [hep- 

ph[. [18 

[8[ A. Czarnecki and K. Melnikov, Nucl. Phys. B544, 520 

(1999); K. G. Chetyrkin, R. Harlander, T. Seidensticker, [19 
and M. Steinhauser, Phys. Rev. D 60, 114015 (1999); 
I. R. Blokland, A. Czarnecki, M. Slusarczyk, and F. Tka- [20 
chov, Phys. Rev. Lett. 93, 062001 (2004); Phys. Rev. D 
71, 054004 (2005); 79, 019901(E) (2009); R. Bonciani [21 
and A. Ferroglia, JHEP 0811, 065 (2008). 

[9] A. Denner and T. Sack, Nucl. Phys. B 358, 46 (1991); [22 
G. Eilam, R. R. Mendel, R. Migneron, and A. Soni, Phys. 
Rev. Lett. 66, 3105 (1991); C.-P. Yuan and T. C. Yuan, [23 
Phys. Rev. D 44, 3603 (1991); T. Kuruma, Z. Phys. C 
57, 551 (1993); S. M. Oliveira, L. Briicher, R. Santos, [24 
and A. Barroso, Phys. Rev. D 64, 017301 (2001). 
[10[ G. Corcella and A. D. Mitov, Nucl. Phys. B 623, 247 [25 
(2002). 

[11[ G. Corcella and F. Mescia, Eur. Phys. J. C 65, 171 

(2010); 68, 687(E) (2010); S. Biswas, K. Melnikov, and [26 
M. Schulze, JHEP 1008, 048 (2010). 

[12| D0 Collaboration, B. Abbott et al., Phys. Rev. D 58, [27 
052001 (1998); 



Phys. Rev. D 63, 
Phys. Rev. 



CDF Collaboration, T. Affolder et al., 
032003 (2001). 

J. Binnewies, B.A. Kniehl, and G. Kramer, 
D 58, 034016 (1998); 

M. Cacciari and M. Greco, Nucl. Phys. B421, 530(1994). 
V. N. Gribov and L. N. Lipatov, Sov. J. Nucl. Phys. 15, 
438 (1972) [Yad. Fiz. 15, 781 (1972)[; G. Altarelli and 
G. Parisi, Nucl. Phys. B126, 298 (1977); Yu. L. Dok- 
shitzer, Sov. Phys. JETP 46, 641 (1977) [Zh. Eksp. Teor. 
Fiz. 73, 1216 (1977)[. 

J. C. CoUins, Phys. Rev. D 66 (1998) 094002. 
J. F. Gunion, H. Haber, G. Kane, and S. Dawson, The 
Higgs Hunter's Guide (Addison- Wesley, Reading, MAA, 
1990), and refrences therein. 

A. Kadeer, J. G. Korner, and M. C. Mauser, Eur. Phys. 
J. C 54, 175 (2008). 

S. M. Moosavi Nejad, "Bottom-hadron production 
through top quark decay" DESY-THESIS-2009-017; 
T. Kinosita, J. Math. Phys. 3, 650 (1962); T. D. Lee and 
M. Nauenberg, Phys. Rev. B 1549 (1964) 133. 

A. Czarnecki and S. Davidson, Phys. Rev. D 47, 3063 
(1993) |arXiv:hep-ph/9208240 . 

J. Liu and Y. P. Yao, Phys. Rev. D 46, 5196 (1992) 
arXiv: hep-ph /9205245 1 . 

B. A. Kniehl, G. Kramer, L Schienbein, and H. Spies- 
berger, Phys. Rev. D 77, 014011 (2008). 

A. Heister et al. (ALEPH Collaboration), Phys. Lett. B 
512, 30 (2001). 

G. Abbiendi et al. (OPAL Collaboration), Eur. Phys. J. 
C 29, 463 (2003). 

K. Abe et al. (SLD Collaboration), Phys. Rev. Lett. 
84, 4300 (2000); Phys. Rev. D 65, 092006 (2002); 66, 
079905(E) (2002). 

Particle Data Group, W. M. Yao et al, J. Phys. G 33, 
(2006) 1. 

J. C. Collins, Phys. Rev. D 58, 094002 (1998). 



